QUARTIC CURVES AND THEIR BITANGENTS 
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Abstract. A smooth quartic curve in the complex projective plane has 36 
inequivalent representations as a symmetric determinant of linear forms and 63 
representations as a sum of three squares. These correspond to Cayley octads 
and Steiner complexes respectively. We present exact algorithms for comput- 
ing these objects from the 28 bitangents. This expresses Vinnikov quartics as 
spectrahedra and positive quartics as Gram matrices. We explore the geometry 
of Gram spectrahedra and we find equations for the variety of Cayley octads. 
Interwoven is an exposition of much of the 19th century theory of plane quartics. 



1. Introduction 

We consider smooth curves in the projective plane defined by ternary quartics 

(1.1) f(x,y,z) = c 400 x' i + c 3W x :i y + c 3 o 1 x 3 z + C22ox 2 y 2 + C2nx 2 yz^ hc 004 ^ 4 , 

whose 15 coefficients Qjk are parameters over the field Q of rational numbers. Our 
goal is to devise exact algorithms for computing the two alternate representations 

(1.2) f(x,y,z) = det(xA + yB + zC), 
where A, B, C are symmetric 4 x 4- matrices, and 

(1.3) f(x,y,z) = qi(x,y,z) 2 + q 2 (x,y,z) 2 + q 3 (x,y,z) 2 } 



where the qi(x,y,z) are quadratic forms. The representation (1.2) is of most 
interest when the real curve Vr(/) consists of two nested ovals. Following Helton- 
Vinnikov [T3] and Henrion |14| . one seeks real symmetric matrices A, B, C whose 



span contains a positive definite matrix. The representation (1.3) is of most 



interest when the real curve Vr(/) is empty. Following Hilbert |16j and Powers- 
Reznick-Scheiderer-Sottile [21], one seeks quadrics qi(x,y,z) with real coefficients. 



We shall explain how to compute all representations (1.2) and (1.3) over C. 

The theory of plane quartic curves is a delightful chapter of 19th century math- 
ematics, with contributions by Aronhold, Cayley, Frobenius, Hesse, Klein, Schot- 
tky, Steiner, Sturm and many others. Textbook references include |H HU |26] . ^ 
started in 1834 with Plucker's result [22j that the complex curve Vc(/) has 28 
bitangents. The linear form i = ax + fly + jz of a bitangent satisfies the identity 

f(x,y,z) = g(x,y,z) 2 + £(x,y,z) ■ h(x,y,z) 

for some quadric g and some cubic h. This translates into a system of polynomial 
equations in (a : : 7), and our algorithms start out by solving these equations. 
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Let K denote the corresponding splitting field, that is, the smallest field exten- 
sion of Q that contains the coefficients a, 0, 7 for all 28 bitangents. The Galois 
group Gal(iT, Q) is very far from being the symmetric group S^s- m fact, if the 
coefficients are general enough, it is the Weyl group of E 7 modulo its center, 



;i-4) 



Gal (if, 1 



W(E 7 )/{±1} = Sp 6 (Z/2Z). 



This group has order 8! -36 = 1451520, and it is not solvable [T2| page 18]. We will 
see a combinatorial representation of this Galois group in Section|3~|(Remark |3.13 ). 
It is based on [TB"t §19] an d Thm. 9]. The connection with Sp 6 (Z/2Z) arises 
from the theory of theta functions jH §5]. For further information see [T2[ §11.4]. 



Naturally, the field extensions needed for (1.2) and (1.3) are much smaller for 



special quartics. As our running example we take the smooth quartic given by 
E(x,y,z) = 25 • (x 4 + y 4 + z 4 ) - 34 



1 2 2 1 2 2 1 2 2\ 

[x y + x z +y z ) 



We call this the Edge quartic. It is one of the curves in the family studied by 



William L. Edge in [7, §14], and it admits a matrix representation (1.2) over 

/ x + 2y 2x + z y-2z\ 

x + 2y y + 2z -2x + z 

2x + z y + 2z x — 2y 

\ y — 2z —2x + z x — 2y 

The sum of three squares representation ( 1.3[ ) is derived from the expression 



1.5) E(x,y,z) 



det 



/ 



;i-6) ( 



x 2 y 2 z 2 xy xz yz) 



25 
-55/2 
-55/2 





21 



-55/2 
25 
25 









-55/2 
25 
25 













21 
-21 







-21 
21 




21 \ 





-84/ 



(x 2 \ 

y 2 

z 2 
xy 
xz 

\y z J 



by factoring the above rank-3 matrix as H T -H where if is a complex 3x6- matrix. 
The real quartic curve Vr(-E') consists of four ovals and is shown in Figure [TJ 



Each of the 28 bitangents of the Edge quartic is defined over Q, but the four 
shown on the right in Figure [TJ are tangent at complex points of the curve. The 
following theorem and Table [TJ summarize the possible shapes of real quartics. 

Theorem 1.7. There are six possible topological types for a smooth quartic curve 
Vr(/) in the real projective plane. They are listed in the first column of Tfo&/e[7J 
Each of these six types corresponds to only one connected component in the com- 
plement of the discriminant A in the 14- dimensional projective space of quartics. 



The classification result in Theorem 1.7 is due to Zeuthen [311. An excellent 



exposition can be found in Salmon's book [26, Chapter VI]. Klein [T7j §5] proved 
that each type is connected in the complement of the discriminant {A = 0}. We 
note that A is a homogeneous polynomial of degree 27 in the 15 coefficients 
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FIGURE 1. The Edge quartic and some of its 28 bitangents 



The real curve 


Cayley octad 


Real bitangents 


Real Steiner complexes 


4 ovals 


8 real points 


28 


63 


3 ovals 


6 real points 


16 


31 


2 non-nested ovals 


4 real points 


8 


15 


1 oval 


2 real points 


4 


7 


2 nested ovals 


real points 


4 


15 


empty curve 


real points 


4 


15 



TABLE 1. The six types of smooth quartics in the real projective plane. 



of /. As a preprocessing step in our algorithms, we use the explicit formula for A 
given in [27, Proposition 6.5] to verify that a given quartic curve Vc(/) is smooth. 

The present paper is organized as follows. In Section 2 we present an algorithm, 
based on Dixon's approach [3J, for computing one determinantal representation 



(1.2). The resulting 4x4-matrices A, B and C specify three quadratic surfaces in 
P 3 whose intersection consists of eight points, known as a Cayley octad. 

In Section 3 we use Cayley octads to compute representatives for all 36 inequiv- 



alent classes of determinantal representations (1.2) of the given quartic /. This 
is accomplished by a combinatorial algorithm developed by Hesse in [T3], which 
realizes the Cremona action [5J on the Cayley octads. The output consists of 36 



symmetric 8x8-matrices (3.4). These have rank 4 and their 28 entries are linear 
forms defining the bitangents. 

In Section 4 we focus on Vinnikov quartics, that is, real quartics consisting of 
two nested ovals. Helton and Vinnikov [13J proved the existence of a representation 
(1.2) over M. We present a symbolic algorithm for computing that representation 
in practice. Our method uses exact arithmetic and writes the convex inner oval 
explicitly as a spectrahedron. This settles a question raised by Henrion [HI §1.2]. 
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In Section 5 we identify sums of three squares with Steiner complexes of bitan- 
gents, and we compute all 63 Gram matrices, i.e. all 6x6-matrices of rank 3 as in 
( 1.6), again using only rational arithmetic over K. This ties in with the results of 
Powers, Reznick, Scheiderer and Sottile in [24J, where it was proved that a smooth 
quartic / has precisely 63 inequivalent representations as a sum of three squares 
(1.3). They strengthened Hilbert's theorem in [16J by showing that precisely eight 
of these 63 are real when / is positive. 

Section [6] is devoted to the boundary and facial structure of the Gram spectra- 
hedron. This is the six-dimensional spectrahedron consisting of all sums of squares 
representations of a fixed positive ternary quartic /. We show that its eight spe- 
cial vertices are connected by 12 edges that form two complete graphs K 4 . We 
also study the structure of the associated semidefinite programming problems. 

Section [7] is devoted to the variety of Cayley octads |5l §LX.3]. We discuss its 
defining equations and its boundary strata, we compute the discriminants of ( 1.2) 
and (1.3), and we end with a classification of nets of real quadrics in P 3 . 

We have implemented most of the algorithms presented in this paper in the 
system SAGEQ Our software and supplementary material on quartic curves and 
Cayley octads can be found at math.berkeley . edu/~cvinzant/quartics .html. 



2. Computing a Symmetric Determinantal Representation 
We now prove, by way of a constructive algorithm, that every smooth quartic 



admits a symmetric determinantal representation (1.2). First we compute the 28 



bitangents, £ = ax + f5y + 72: . Working on the affine chart {7 = 1}, we equate 

f(x, y, -ax - (3y) = (k x 2 + n 4 xy + K 2 y 2 ) 2 , 

eliminate Ko,Ki,/«2, and solve the resulting system for the unknowns a and /3. 
This constructs the splitting field K for the given / as a finite extension of Q. All 
further computations in this section are performed via rational arithmetic in K. 

Next consider any one of the ( 2 3 8 ) = 3276 triples of bitangents. Multiply their 
defining linear forms. The resulting polynomial t>oo — ^1^3 is a contact cubic 
for Vc(/)j which means that the ideal (vqq, f) in K[x,y,z] defines six points in 
P 2 each of multiplicity 2. Six points that span three lines in P 2 impose inde- 
pendent conditions on cubics, so the space of cubics in the radical of (vqq, f) is 
4-dimensional over K. We extend {v 00} to a basis {vqq, ^01 > v 02i ^03} of that space. 

Max Noether's Fundamental Theorem [9, § 5.5] can be applied to the cubic v 00 
and the quartic /. It implies that a homogeneous polynomial lies in (vqq, /) if it 
vanishes to order two at each of the six points of Vc({voo, /))■ The latter property 
holds for the sextic forms v Qi Voj. Hence foi^oj lies in (vqq, f) for 1 < i < j < 3. 
Using the Extended Buchberger Algorithm, we can compute cubics Vij such that 

(2.1) v 0i v 0j - v 00 Vij E (/). 



1 www . sagemath . org 
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We now form a symmetric 4x4-matrix V whose entries are cubics in K[x, y, z\: 

fvoo voi V02 V 03 \ 

y VOI «11 Vv2 «13 

V02 V 12 V 22 V 23 

\V03 V 13 V 23 V33J 

The following result is due to Dixon [3J, and it almost solves our problem. 
Proposition 2.2. Each entry of the adjoint V ad i is a linear form times f 2 , and 
det(f~ 2 • V ad j) = 7 • f[x, y, z) for some constant 7 G K . 



Hence, if det(V) 7^ then f 2 ■ V 3 *^ gives a linear matrix representation (1.2) 



Proof. Since vqq G" (/), the condition (2.1) implies that, over the quotient ring 



K[x, y, z]/{f), the matrix V has rank 1. Hence, in the polynomial ring K[x, y, z], 
the cubic / divides all 2 x 2 minors of V. This implies that f 2 divides all 3 x 3 
minors of V, and f 3 divides det(V). As the entries of V ad i have degree 9, it follows 
that = f 2 -W, where W is a symmetric matrix whose entries are linear forms. 
Similarly, as det(V) has degree 12, we have det(V) = 8f 3 for some 5 G K, and 
5 7^ unless det(V) is identically zero. Let I4 denote the identity matrix. Then 

5f-h = det(V)-h = V-V^ = f 2 -V-W. 

Dividing by f 2 and taking determinants yields 

S A f = det(V) • det(W) = Sf 3 ■ det{W). 

This implies the desired identity det(W) = 5 3 f. □ 

We now identify the conditions to ensure that det(V) is not the zero polynomial. 

Theorem 2.3. The determinant of V vanishes if and only if the six points of 
Vc(/j ^1^3), at which the bitangents £1,^2,^3 touch the quartic curve Vc(f), He 
on a conic in P 2 . This happens for precisely 1260 of the 3276 triples of bitangents. 

Proof. Dixon [3J proves the first assertion. The census of triples appears in the 
table on page 233 in Salmon's book [26, §262]. It is best understood via the Cayley 
octads in Section 3. For further information see Dolgachev's notes |U §6.1]. □ 

Remark 2.4. Let £1, £2, £3 be any three bitangents of Vc(/)- If the six intersection 
points with Vc(f) he on a conic, the triple {£1,^2,^3} is called syzygetic, otherwise 
azygetic. A smooth quartic / has 1260 syzygetic and 2016 azygetic triples of 
bitangents. Similarly, a quadruple {£1, £2, £ 3 , £4} of bitangents is called syzygetic if 
its eight contact points lie on a conic and azygetic if they do not. Every syzygetic 
triple £i,£2,£3 determines a fourth bitangent £4 with which it forms a syzygetic 
quadruple. Indeed, if the contact points of £1,^2,^3 he on a conic with defining 
polynomial q, then q 2 lies in the ideal (j,t\lilz)i so that q 2 = 7/ + l\l<il 3 l^ and 
the other two points in Vc(/, q) must be the contact points of the bitangent £4. 
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Algorithm 2.5. Given a smooth ternary quartic / G Q[x,y, z], we compute the 
splitting field K over which the 28 bitangents of Vc(/) are defined. We pick a 
random triple of bitangents and construct the matrix V via the above method. If 
det(y) 7^ 0, we compute the adjoint of V and divide by f 2 , obtaining the desired 
determinantal representation of / over K. If det(V) = 0, we pick a different triple 
of bitangents. On each iteration, the probability for det(V r ) 7^ is |||| = ^. 

Example 2.6. The diagram on the left of Figure [T] shows an azygetic triple of 
bitangents to the Edge quartic. Here, the six points of tangency do not lie on a 
conic. The representation of the Edge quartic in ( |1.5[ ) is produced by Algorithm 2.5 



starting from the contact cubic t>oo = 2(y + 2z)(— 2x + z)(x — 2y). □ 
3. CAYLEY OCTADS AND THE CREMONA ACTION 



Algorithm 2.5 outputs a matrix M = xA+yB+zC where A, B, C are symmetric 



4x4- matrices with entries in the subfield K of C over which all 28 bitangents of 



Vc(/) are defined. Given one such representation (1.2) of the quartic /, we shall 



construct a representative from each of the 35 other equivalence classes. Two 



representations (1.2) are considered equivalent if they are in the same orbit under 
the action of GL 4 (C) by conjugation M 1— > U T MU. We shall present an algorithm 
for the following result. It performs rational arithmetic over the splitting field K 
of the 28 bitangents, and it constructs one representative for each of the 36 orbits. 

Theorem 3.1 (Hesse [15J). Every smooth quartic curve f has exactly 36 equiva- 



lence classes of linear symmetric determinantal representations (1.2). 

Our algorithms begins by intersecting the three quadric surfaces seen in M: 

(3.2) uAu T = uBu T = uCu T = where u = (u : u x : u 2 : u 3 ) G P 3 (C). 

These equations have eight solutions 0\, . . . ,0%. This is the Cayley octad of M. 
In general, a Cayley octad is the complete intersection of three quadrics in P 3 (C). 

The next proposition gives a bijection between the 28 bitangents of Vc(/) and 
the lines OiOj for 1 < i < j < 8. The combinatorial structure of this configuration 
of 28 lines in P 3 plays an important role for our algorithms. 

Proposition 3.3. Let Oi, ...,0s be the Cayley octad defined above. Then the 28 
linear forms O.MOj G C[x, y,z] are the equations of the bitangents ofVc(f)- 

Proof. Fix i j. After a change of basis on C 4 given by a matrix U G GL^C) and 
replacing M by U T MU, we may assume that Oi = (1, 0, 0, 0) and Oj = (0, 1, 0, 0). 
The linear form b^ = OiMOj now appears in the matrix: 



M 



bij 
hj 


M' 


(M'f 


* 



Expanding det(M) and sorting for terms containing b^ shows that /=det(M) is 
congruent to det(M') 2 modulo (bij). This means that b^ is a bitangent. □ 
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Let O be the 8 x 4-matrix with rows given by the Cayley octad. The symmetric 
8x8-matrix OMO T has rank 4, and we call it the bitangent matrix of M. By the 
definition of O, the bitangent matrix has zeros on the diagonal, and, by Proposi- 



tion 3.3, its 28 off-diagonal entries are precisely the equations of the bitangents: 



(3.4) 



OMO 1 



(o 


bn 


bl3 


&14 


hr, 


&16 


bn 




b\2 





b23 


&24 


b25 


&26 


b27 


^28 


bn 


b23 





&34 


b35 


&36 


b37 


^38 


bu 


b24 


&34 





&45 


&46 


647 


^48 


h5 


b25 


&35 


&45 





&56 


&57 


^58 


he 


b26 


&36 


bi6 


^56 





b&7 


^68 


bn 


b27 


&37 


647 


b57 


&67 





^78 


\bis 


b28 


&38 


bis 


b58 


bes 


&78 


0/ 


at the octad 


1; ■ ■ 


,o s 


consists of K- 



of P 3 : To see this, let if' be the field of definition of the octad over K. Then any 
element a of Gal (if : K) acts on the octad by permutation, and thus permutes 
the indices of the bitangents, fey. On the other hand, as all bitangents are defined 
over K, a must fix fey (up to a constant factor). Thus the permutation induced 
by a on the octad must be the identity and Gal(K' : K) is the trivial group. 

Example 3.6. The symmetric matrix M in ( 1.5[ ) determines the Cayley octad 

/l -1 1 1 3 \ 



J 



10 3 -11 1 
10 1 3 1-1 
\0 1 -1 -1 3 -1/ 

All the 28 bitangents of E(x, y, z) are revealed in the bitangent matrix OMO T = 

t + y+z \ 



by — 5z — 3x 
24x + 12y 
12x + 24z 
242/ - 12z 




/ 


x + 2y 


2x + z 


y-2z 


5x+5y+3z 


5x — 3y+5z 


3x + by — 5z 





y + 2z 


-2x + z 


x — y+z 


3x+5y+5z 


— 5x+3y+5z 


2x + z 


y + 22 





x-2y 


— 3x+5z+5y 


x — z+y 


5x + 3z — 5y 


y-2z 


-2x + z 


!C-2j/ 





—3y+5z — 5x 


—5y — 3z + 5x 


—y—z— x 


5x+5y+3z 


x-y *-z 


— 3x+5z+5j/ 


—3y+5z — 5x 





24y + 12z 


-I2x + 24z 


5x — 3y+5z 


3x+5y+5z 


x — z+y 


— 5y — 3z+5x 


24y + 12z 





24x - 12y 


3x+5y — 5z 
\ —x + y+z 


— 5x+3y+5z 


5x+3z — 5y 


-y-z-x 


-12a + 24z 


24a; - 12y 





5x+5y—3z 


5x+5z+3y 


5y — 5z — 3x 


24x + 12y 


I2x + 24z 


24y - Viz 



Each principal 4x4-minors of this matrix is a multiple of E(x, y, z), as in (3.7). □ 



Each principal 3x3-minor of the bitangent matrix ( 3.4[ ) is a contact cubic 
2bijbikbjk of Vc(jf) and can serve as the starting point for the procedure in Section 
2. Hence, each principal 4x4- minor My-jw of (3.4) represents the same quartic: 



, „ _s det(Mijki) = a non- zero scalar multiple of f(x, y, z) 

= b\b\ x + b\ k b\ + b\b) k - 2(b ij b ik b j ib M + b i:j b a b jk b k i + b^bub^bji) 



However, all these = 70 realizations of (1.2) lie in the same equivalence class. 

In what follows, we present a simple recipe due to Hesse [TJ] for finding 35 
alternate bitangent matrices, each of which lies in a different GL 4 (C)-orbit. This 
furnishes all 36 inequivalent determinantal representations promised in Theorem 



3.1 We begin with a remark that explains the number 1260 in Theorem 2.3 
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Remark 3.8. We can use the combinatorics of the Cayley octad to classify 
syzygetic collections of bitangents. There are 56 triples A of the form {bij, b^, bj^}. 



Any such triple is azygetic, by the if-direction in Theorem |2.3[ because the cubic 

Ev- 



bijbikbjk appears on the diagonal of the adjoint of the invertible matrix M, 



ijkl- 



ery product of an azygetic triple of bitangents appears as a 3 x 3 minor of exactly 
one of the 36 inequivalent bitangent matrices, giving 36-56 = 2016 azygetic triples 
of bitangents and ( 2 3 8 ) — 2016 = 1260 syzygetic triples. 

A quadruple of bitangents of type □ is of the form {b^, bjk, b^i, bu}. Any such 



quadruple is syzygetic. Indeed, equation (3.7) implies f-\-^{hjbjkhihi) = (hjbki — 



hkbji + bubjk) , and this reveals a conic containing the eight points of contact. 



Consider the following matrix which is gotten by permuting the entries of M, 



ijkl ■ 



M' ijkl 



/0 

hi 
bji 
\bjk 



bu 


bu 
bik 



bji 
hi 


bij 



bjk\ 
bik 
bij 
0/ 



This procedure does not change the determinant: det(M^- fcZ ) = det(Mijki) = f . 
This gives us 70 linear determinantal representations (1.2) of the quartic /, one 
for each quadruple I = {i,j,k,l} C {1, . . . , 8}. These are equivalent in pairs: 

Theorem 3.9. If I ^ J are quadruples in {1, ... , 8}, then the symmetric matrices 
M\ and M'j are in the same GL^(C)-orbit if and only if I and J are disjoint. None 
of these orbits contains the original matrix M = xA + yB + zC . 

Proof. Fix / = {1, 2, 3, 4} and note the following identity in K[x, y, z, uq, Ui, U2, U3]: 



MA7 



Ui 

u 2 

w 



bi2 
bu 
\bu 



bu 


&23 
&24 



6l3 



b 3 4 



bu\ 

0/ 



f u o\ 

w 



/V\7 



u u 1 u 2 u 3 



Ui 1 

u 2 1 



&34 
&24 
\&23 



&34 


b u 
bi3 



b 2 4 
bu 


bi2 



b23\ 

h 3 

0/ 



u 2 1 



This shows that the Cayley octad of M[ 234 is obtained from the Cayley octad of 
M1234 by applying the Cremona transformation at 0\, 2 , 03,0/t. Equivalently, 
observe that the standard basis vectors of Q 4 are the first four points in the Cayley 



octads of both M1234 and M{ 234 , and if O 



[Oti 



[OLi 



: Si) 

&T 1 



for 



5,6,7,8 
for 



belong to the Cayley octad of M1234, then 0\ 
i = 5, 6, 7, 8 belong to the Cayley octad O' of M( 234 . 

Thus the transformation from Mijki to M[- kl corresponds to the Cremona action 
cr3 g on Cayley octads, as described on page 107 in the book of Dolgachev and 
Ortland [5j. Each Cremona transformation changes the projective equivalence 
class of the Cayley octad, and altogether we recover the 36 distinct classes. That 
M'j is equivalent to M'j when I and J are disjoint can be explained by the following 
result due to Coble [2]. See [51 §111.3] for a derivation in modern terms. □ 

Theorem 3.10. Let O be an unlabeled configuration of eight points in linearly 
general position in P 3 . Then O is a Cayley octad (i.e. the intersection of three 
quadrics) if and only if O is self- associated (i.e. fixed under Gale duality; cf. 
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The Cremona action on Cayley octads was known classically as the bifid substi- 
tution, a term coined by Arthur Cayley himself. We can regard this as a combi- 
natorial rule that permutes and scales the 28 entries of the 8x8 bitangent matrix: 

Corollary 3.11. The entries of the two bitangent matrices OM 1234 T = (6y) and 
0'M' l2M 0' T = (b'^) are related by non-zero scalars in the field K as follows: 

f b kl if {i,j,k,l} = {1,2,3,4}, 
The linear form b'^ is a scalar multiple of < fry if |{«, j} D {1,2,3,4}| = 1, 

( b kl if {i,j, k, 1} = {5, 6, 7, 8}. 

Proof. The first case is the definition of M{ 234 . For the second case we note that 

m o\ & 15 = OiM l2M Ol = (3 5 b 12 + 75^13 + #5&14 

[6 - iA) and b' 15 = 0[M[ 2U 0' 5 T = fc l b M + j^b 24 + S^b 23 , 



by Proposition 3.3 The identity O5M1234O5 = 0, when combined with (3.12), 
translates into 05615 + (3^5^'^ = 0, and hence b' 15 = —a 5 (3^ 1 'j^ 1 5^ 1 bi5. F° r the 
last case we consider any pair {i,j} C {5,6,7,8}. We know that b'^ = vb k i, for 
some v G K* and {k, 1} C {5, 6, 7, 8}, by the previous two cases. We must exclude 
the possibility {k,l} R {i,j} 7^ 0. After relabeling this would mean b' 56 = vb§§ or 
b' 56 = ub 57 . If b' 56 = z/6 56 then the lines {b' 12 , b' 25 , b' 56 , b' w } and {6 34 , b 25 , b 5e , b ie } 
coincide. This is impossible because the left quadruple is syzygetic while the right 



quadruple is not, by Remark 3.8 Likewise, 6' 56 = z/657 would imply that the 



azygetic triple {6' 15 , b' m , b' 16 } corresponds to the syzygetic triple {615, 6 57 , &i 6 }. □ 

Remark 3.13. The 35 bifid substitutions of the Cayley octad are indexed by 
partitions of [8] = {1, 2, . . . , 8} into pairs of 4-sets. They are discussed in modern 
language in [5, Prop. 4, page 172]. Each bifid substitution determines a permu- 
tation of the set (®) = {{i,j} ■ l<i<7<8}. For instance, the bifid partition 



1234|5678 determines the permutation in Corollary 3.11 Hesse [T51 page 318] 
wrote these 35 permutations of 0®) explicitly in a table of format 35x28. Hesse's 



remarkable table is a combinatorial realization of the Galois group (1.4). Namely 



W(E 7 )/{±1} is the subgroup of column permutations that fixes the rows. 

We conclude this section with a remark on the real case. Suppose that / is given 



by a real symmetric determinantal representation (1.2), i.e. / = det(M) where 
M = xA + yB + zC and A, B, C are real symmetric 4 x 4- matrices. By [301 §0], 
such a representation exists for every smooth real quartic /. Then the quadrics 
uAu T ,uBu T ,uCu T G K[uo,ui,u 2 ,us} 2 defining the Cayley octad are real, so that 
the points 0\, ... ,Os are either real or come in conjugate pairs. 

Corollary 3.14. Let M = xA+yB+zC be a real symmetric matrix representation 
of f with Cayley octad Oi, ... ,Os- Then the bitangent OfMOj is defined over M. 
if and only if Oi and Oj are either real or form a conjugate pair, Oi = Oj . 

From the possible numbers of real octad points we can infer the numbers of 
real bitangents stated in Table [TJ If 2k of the eight points are real, then there are 
4 — k complex conjugate pairs, giving ( 2 2 fc ) + 4 — k = 2k 2 — 2k + 4 real bitangents. 
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FIGURE 2. The Vinnikov quartic in Example 4.1 



4. Spectrahedral Representations of Vinnikov Quartics 

The symmetric determinantal representations / = det(M) of a ternary quartic 
/ G Q[x,y,z] are grouped into 36 orbits under the action of GL^C) given by 
M i — y T T MT. The algorithms in Sections [2] and [3] construct representatives for 



all 36 orbits. If we represent each orbit by its 8x8-bitangent matrix (3.4), then this 
serves as a classifier for the 36 orbits. Suppose we are given any other symmetric 
linear matrix representation M = xA + yB + zC of the same quartic /, and our 
task is to identify in which of the 36 orbits it lies. We do this by computing the 
Cayley octad O of M and the resulting bitangent matrix OMO T . That 8x8- 
matrix can be located in our list of 36 bitangent matrices by comparing principal 
minors of size 3x3. These minors are products of azygetic triples of bitangents, 
and they uniquely identify the orbit since there are 2016 = 36-56 azygetic triples. 

We now address the problem of finding matrices A, B and C whose entries are 
real numbers. Theorem 11.71 shows that this is not a trivial matter because none 



of the 36 bitangent matrices in (3.4) has only real entries, unless the curve Vr(/) 



consists of four ovals (as in Figure pj. We discuss the case when the curve is a 
Vinnikov quartic, which means that Vr(/) consists of two nested ovals. 
As shown in [13], the region bounded by the inner oval corresponds exactly to 

{(x,y, z) E 1R 3 : xA + yB + zC is positive definite}, 

a convex cone. This means that the inner oval is a spectrahedron. The study of such 
spectrahedral representations is of considerable interest in convex optimization. 
Recent work by Henrion [2] underscores the difficulty of this problem for curves 
of genus g > 2, and in the last two paragraphs of [T4"] §1.2], he asks for the 
development of a practical implementation. This section constitutes a definitive 
computer algebra solution to Henrion's problem for smooth quartic curves. 

Example 4.1. The following smooth quartic is a Vinnikov curve: 
f{x,y,z) = 2x 4 + y 4 + z 4 - 3x 2 y 2 - 3xV + y 2 z 2 . 
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Running the algorithm in Section 2, we find that the coefficients of the 28 bitan- 
gents are expressed in radicals over Q. However, only four of the bitangents are 
real. Using Theorem |4. 3 below, we conclude that there exists a real matrix repre- 



sentation 1.2 with entries expressed in radicals over Q. One such representation is 

(ux + y az bz \ 

ux — y cz dz 

az cz x + y 

\ bz dz x — yj 



(4.2) f(x,y,z) = det 



with 



a = -0.57464203209296160548032752478263071485849363449367..., 

b = 1.03492595196395554058118944258225904539129257996969..., 

c = 0.69970597091301262923557093892256027951096114611925..., 

d = 0.4800486503802432010856027835498806214572648351951..., 

u = V2 = 1.4142135623730950488016887242096980785696718.... 

The expression in radicals is given by the following maximal ideal in Q[a, b, c, d, u]: 

(u 2 - 2 , 256d 8 - 384d 6 u+256d 6 -384d 4 w+672d 4 -336d 2 u+448d 2 -84u+121, 
23c + 7584c/ 7 M+10688rf 7 -5872c/ 5 M-8384c/ 5 +1806c/ 3 M+2452rf 3 -181rfM-307c/, 
236 + 5760rf 7 u+8192rf 7 -4688rf 5 u-6512d 5 +1452d 3 u+2200rf 3 -212rfM-232(i, 
23a - 1440rf 7 u-2048d 7 +1632rf 5 u+2272d 5 -570d 3 M-872rf 3 +99du+81d). 

A picture of the curve Vr(/) in the affine plane {x = 1} is shown in Figure [2j □ 



The objective of this section is to establish the following algorithmic result: 

Theorem 4.3. Let f £ Q[x,y,z] be a quartic whose curve Vc(/) is smooth. 
Suppose /(x, 0,0) = x 4 and f(x,y,0) is squarefree, and let K be the splitting field 
for its 28 bitangents. Then we can compute a determinantal representation 

(4.4) f(x,y,z) = det(xl + yD + zR) 

where I is the identity matrix, D is a diagonal matrix, R is a symmetric matrix, 
and the entries of D and R are expressed in radicals over K . Moreover, there 
exist such matrices D and R with real entries if and only if Vr(/) is a Vinnikov 
curve containing the point (1:0:0) inside the inner oval. 

The hypotheses in Theorem |4.3| impose no loss of generality. Any smooth 
quartic will satisfy them after a linear change of coordinates (x : y : z) in P 2 . 

Proof. Using the method in Section 2, we find a first representation f(x,y,z) = 
det(xA + yB + zC) over the field K. However, the resulting matrices A,B,C 
might have non-real entries. The matrix A is invertible because we have assumed 
det(xA) = f(x, 0,0) = x 4 , which implies det (A) = 1. 

The binary form f(x,y,0) = det(xA + yB) is squarefree. That assumption 
guarantees that the 4x4-matrix A~ 1 B has four distinct complex eigenvalues. Since 
its entries are in K, its four eigenvalues lie in a radical extension field L over K. 
By choosing a suitable basis of eigenvectors, we find a matrix U £ GL^L) such 
that U~ 1 A~ 1 BU is a diagonal matrix D 1 = diag(Ai, A 2 , A 3 , A 4 ) over the field L. 
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We claim that D 2 = U T AU and D 3 = U T BU are diagonal matrices. For each 
column Ui of U we have A~ 1 Bui = AjU, , so Bui = AjAu,. For 1 < i < j < 4 this 
implies ujBui = XiujAu^ and, by switching indices, we get ujBuj = XjujAuj. 
Since B is symmetric, the difference of the last two expressions is zero, and we 
conclude (Aj — Xj) ■ ujAuj = 0. By assumption, we have Aj 7^ Aj and therefore 
ujAuj = and ujBuj = 0. This means that D 2 and D 3 are diagonal. 

Let D 4 be the diagonal matrix whose entries are the reciprocals of the square 
roots of the entries of D 2 - These entries are also expressed in radicals over K. 
Then D 4 D 2 D 4 = I is the identity matrix, D 4 D 3 D 4 = D is also diagonal, and 



D 4 U T MUD 4 = xI + yD + zR 



is the real symmetric matrix representation required in (4.4). 

In order for the entries of D and R to be real numbers, it is necessary (by |13j ) 
that Vr(/) be a Vinnikov curve. We now assume that this is the case. The exis- 



tence of a real representation (4.4) is due to Vinnikov |3"Ut §0]. A transcendental 
formula for the matrix entries of D and R in terms of theta functions is pre- 
sented in equations (4.2) and (4.3) of [13, §4]. We need to show how our algebraic 
construction above can be used to compute Vinnikov's matrices D and R. 

Given a quartic / G Q[x, y, z] with leading term x A , the identity (4.4) translates 



into a system of 14 polynomial equations in 14 unknowns, namely the four entries 
of D and the ten entries of R. For an illustration of how to solve them see Example 



4.6 We claim that these equations have at most 24-8-36 = 6912 complex solutions 



and all solutions are expressed in radicals over K. Indeed, there are 36 conjugation 



orbits, and per orbit we have the freedom to transform (4.4) by a matrix T such 
that T T T = I and T T DT is diagonal. Since the entries entries of D are distinct, 
these constraints imply that T is a permutation matrix times a diagonal matrix 
with entries ±1. There are 24 • 16 possible choices for T, but T and — T yield the 
same triple (J, D, R), so the number of solutions per orbit is 24 • 8. 



We conclude that, for each of the 36 orbits, either all representations (4.4) are 
real or none of them is. Hence, by applying this method to all 36 inequivalent 
symmetric linear determinantal representations constructed in Section 3, we are 
guaranteed to find Vinnikov's real matrices D and R. See also |2T1 Section 2] for 
additional examples and a more detailed discussion. □ 

The above argument for the simultaneous diagonalizability of A and B is taken 
from Greub's linear algebra text book [11J. We could also handle the exceptional 
case when A _1 B does not have four distinct eigenvalues. Even in that case there 
exists a matrix U in radicals over K such that U T AU and U T BU are diagonal, 
but the construction of U is more difficult. The details are found in [Til §IX.3]. 

Corollary 4.5. Every smooth Vinnikov curve has a real determinantal represen- 



tation (1.2) in radicals over the splitting field K of its 28 bitangents. 



We close with the remark that the representation (4.4) generally does not exist 
over the field K itself but the passage to a radical extension field is necessary. 
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Example 4.6. All 6912 matrix representations xl + yD + zR of the Edge quartic 
E(x, y, z) = 25-(x A +y 4 +z 4: ) — 34-(x 2 y 2 +x 2 z 2 +y 2 z 2 ) are non-real and have degree 4 
over Q. The entries of D are the four complex zeros of the irreducible polynomial 
x 2 + 1. After fixing D, we have 192 choices for R, namely, selecting one of 



x 



25 



the 36 orbits fixes R up to conjugation by diag(±l, ±1, ±1, ±1). For the orbit of 



the matrix xA+yB+zC in (1.5), our algorithm gives the representation 

/-v / 21/5 - 2t/5 




D 



2i/5 



V 













R 




3/7 + i) 




^27/35 



5. Sums of Three Squares and Steiner Complexes 




□ 



Our next goal is to write the given quartic / as the sum of three squares of 



quadrics. Such representations (1.3) are classified by Gram matrices of rank 3. A 



Gram matrix for / is a symmetric 6x6 matrix G with entries in C such that 



./ 



• G ■ v where 



/ 2 2 2 \j 

[x ,y ,z ,xy,xz,yz) 



v 

We can write G = H T ■ H, where H is an r x 6-matrix and r = rank(G). Then 
the factorization / = (Hv) T ■ (Hv) expresses / as the sum of r squares. 

It can be shown that no Gram matrix with r < 2 exists when / is smooth, and 
there are infinitely many for r > 4. For r = 3 their number is 63 by Theorem |5.1 



Gram matrices classify the representations (|1.3|): two distinct representations 
/ 



2 , 2 , 2 
?1 + Q2 + % 



2 , 2 , 2 
Pl+P2+ P3 



correspond to the same Gram matrix G of rank 3 if and only if there exists an 
orthogonal matrix T G 03(C) such that T • (pi,P2,P3) T = (<?i, Q2, <?3) T - The 
objective of this section is to present an algorithmic proof for the following result. 

Theorem 5.1. Let f G Q[x,y,z] be a smooth quartic and K the splitting field for 
its 28 bitangents. Then f has precisely 63 Gram matrices of rank 3, all of which 
we compute using rational arithmetic over the field K . 

The fact that / has 63 Gram matrices of rank 3 is a known result due to Coble 
(21 Ch. 1, §14]; see also f2M Prop. 2.1]. Our contribution is a new proof that 
yields a i^-rational algorithm for computing all rank-3 Gram matrices. Instead of 
appealing to the Jacobian threefold of /, as in [24J, we shall identify the 63 Gram 
matrices with the 63 Steiner complexes of bitangents (see [2"6l §VI] and [4, §6]). 

We begin by constructing a representation / = q\ + q\ + q\ from any pair of 
bitangents. Let £,£' be distinct bitangents of /, and let p G C[x,y,z] 2 be a non- 
singular quadric passing through the four contact points of IE' with /. By Max 
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Noether's Fundamental Theorem [9] § 5.5], the ideal (££', /) contains p 2 , thus 
(5.2) / = ££'u-p 2 , 

for some quadric u G C[x,y, z] 2 , after rescaling p by a constant. Over C, the 



identity (5.2) translates directly into one of the form: 



(5.3) / = g«< + ^)\(I«<-I u ) 2 + (^. 

Remark 5.4. Just as systems of contact cubics to Vc(/) were behind the formula 
, systems of contact conies to Vc(/) are responsible for the representations 
. The simplest choice of a contact conic is a product of two bitangents. 



1.2 



1.3 



In (5.3) we wrote / as a sum of three squares over C. There are ( 2 2 8 ) = 378 pairs 



{£,£'} of bitangents. We will see Theorem 5.10 that each pair forms a syzygetic 
quadruple with 5 other pairs. This yields 378/6 = 63 equivalence classes. More 
importantly, there is a combinatorial rule for determining these 63 classes from a 
Cayley octad. This allows us to compute the 63 Gram matrices over K. 



Equation (]5.2|) can also be read as a quadratic determinantal representation 

Qo Qi 



(5.5) / = det 

\Ql 12 

with go = ££', Qi — P, and g 2 = u. This expression gives rise to the quadratic 
system of contact conies {\lq Q + 2\ Q \ l q 1 + \\ql : A G P 1 (C)}. The implicitization 
of this quadratic system is a quadratic form on span{go, Qi, 92}- With respect to 
the basis (go, Qi, Q2), it is represented by a symmetric 3x3 matrix C. Namely, 




C — I —1 and its inverse is C 



-1 




The formula (J5T5J) shows that / = g g 2 - q\ — (go, gi, g 2 ) • C 1 ■ (g , gi, g 2 ) T . We 
now extend g , gi, g 2 to a basis g = (g , gi, g 2 , g3, Qi, Qb) °f y, z} 2 . Let T denote 
the matrix that takes the monomial basis v = (x 2 ,y 2 , z 2 ,xy,xz,yz) to g. If G is 
the 6x6 matrix with C~ l in the top left block and zeros elsewhere, then 

(5.6) / = (go,gi,g 2 )-C- 1 -(go,g 1 ,g 2 ) T = v T ■ T T ■ G ■ T ■ v. 

Thus, G = T T GT is a rank-3 Gram matrix of /. This construction is completely 
reversible, showing that every rank-3 Gram matrix of / is obtained in this way. 



The key player in the formula (5.6) is the quadratic form given by C. From this, 
one easily gets the Gram matrix G. We shall explain how to find G geometrically 
from the pair of bitangents £,£'. The following result is taken from Salmon [26J: 

Proposition 5.7. Let f = det(Q) where Q is a symmetric 2 x 2-matrix with 



entries in C[x,y,z]2 as in (5.5). Then Q defines a quadratic system of contact 



conies \ T Q\, A G P 1 (C) ; that contains exactly six products of two bitangents. 
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Sketch of Proof. To see that X T QX is a contact conic, note that for any A, fi G C 2 , 
(5.8) (\ T Q\)(fi T Qfi) - (\ T Qfif = ^iVkMQijQki ~ QikQji). 



The expression QjjQ 



kl 



i,j,k,l 

QikQji is a multiple of det(Q) 



/, and hence so is the 



left hand side of (5.8). This shows that X T QX is a contact conic of Vc(/)- The 
set of singular conies is a cubic hypersurface in C[x,y, z\%. As X T QX is quadratic 
in A, we see that there are six points A G P X (C) for which X T QX is the product of 
two linear forms. These are bitangents of / and therefore ^-rational. □ 

Remark 5.9. If the Gram matrix G is real, then it is positive (or negative) 
semidefinite if and only if the quadratic system Q = {X T QX \ X G P 1 (C)} does 
not contain any real conies. For if G is real, we may take a real basis (q' , q[, q' 2 ) 
of span{go) Qi, Q2} = ker(G)- 1 in C[x, y,z\i- If Q does not contain any real conies, 
then the matrix C' representing Q with respect to the basis (q' Q , q[, q' 2 ) is definite. 
Using C' instead of C in the above construction, we conclude that C'~ x is definite 
and hence G is semidefinite. The converse follows by reversing the argument. 

We now come to Steiner complexes, the second topic in the section title. 

Theorem 5.10. Let S = {{£\, £[},...,{£$, £' 6 }} be six pairs of bitangents of a 
smooth quartic f G Q[x,y, z\. Then the following three conditions are equivalent: 

(1) The reducible quadrics £\£\ , . . . , £e£' e lie in a system of contact conies X T QX, 
X G P X (C), for Q a quadratic determinantal representation (5.5) of f . 

(2) For each i ^ j , the eight contact points Vc{£i£'i£j£'j) H Vc(/) lie on a conic. 

(3) With indices as in the bitangent matrix (3.4 ) for a Cayley octad, either 

<S = {{bik,bjk} I {i, j} = I and k G I c } for a 2-set I C {1, . . . ,8} ; 

or S = {{bij.bki} I {i,j, k,l} = I or {i,j, k,l} = I c ) for a A-set I C {1, ... ,8}. 

Proof. This is a classical result due to Otto Hesse [T5] . The proof can also be found 
in the books of Salmon [26J and Miller-Blichfeldt-Dickson UB §185-186]. □ 

A Steiner complex is a sextuple S of pairs of bitangents satisfying the conditions 



of Theorem 5.10 A pair of bitangents in S is either of the form {6^, bjk} (referred 
to as type \/) or of the form {6^-, b^} (type ||). The first type of Steiner complex 



in Theorem |5.10 (3) contains pairs of bitangents of type \f and the second type 
contains pairs of type ||. There are (T\ = 28 Steiner complexes of type \f and 
(*)/2 = 35 Steiner complexes of bitangents of type ||. The two types of Steiner 
complexes are easy to remember by the following combinatorial pictures: 






Type V Type 1 1 

This combinatorial encoding of Steiner complexes enables us to derive the last 
column in Table [l] in the Introduction. We represent the quartic as (1.3) with 
A, B, C real, as in [3D]. The corresponding Cayley octad {Ox, . . . , Og} is invari- 
ant under complex conjugation. Let n be the permutation in S% that represents 
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complex conjugation, meaning Oj = O^ny Then complex conjugation on the 63 
Steiner complexes is given by the action of tt on their labels. For instance, when 
all Oi are real, as in the first row of Table [TJ then tt is the identity. For the 
other rows we can relabel so that tt = (12), tt = (12)(34), tt = (12)(34)(56) and 
7i = (12) (34) (56) (78). We say that a Steiner complex S is real if its labels are fixed 
under tt. For example, if S is the Steiner complex {{&13, &23jv • • , {bis, &2s}} of type 
V as above, then S is real if and only if tt fixes {1,2}. Similarly, if S is the Steiner 
complex {{612, 634}, {bi3,b 2 i}, {&i 4 ,&23}, {b 5 6,b 78 }, {657 As}, {b 5 8,b 67 }} of type 
||, then S is real if and only if tt fixes the partition {{1, 2, 3, 4}, {5, 6, 7, 8}}. For 
instance, for the empty curve, in the last row Table [TJ one can check that exactly 
15 Steiner complexes are fixed by tt = (12)(34)(56)(78), as listed in Section [6j 

We now sum up what we have achieved in this section, namely, a recipe for 
constructing the 63 Gram matrices from the 28 + 35 Steiner complexes \f and 1 1 . 



Proof and Algorithm for Theorem 5.1, We take as input a smooth ternary quar- 
tic / G Q[x,y,z] and any of the 63 Steiner complexes {{£i,£'i}, ■ ■ ■ ,{£e,£' 6 }} of 
bitangents of Vc(f). From this we can compute a rank-3 Gram matrix for / as 
follows. The six contact conies £i£\ span a 3-dimensional subspace of K[x,y, z] 2 , 



by Theorem 5.10 (1), of which {£i£i, ti^ii ^3^3} * s a basis. The six vectors £i£\ lie 
on a conic in that subspace, and we compute the symmetric 3x3-matrix C repre- 
senting this conic in the chosen basis. We then extend its inverse C~ l by zeroes 
to a 6x6 matrix G and fix an arbitrary basis {54, q$, q$} of spaa{^i^, £2^2, £3^3} 
in K[x,y, z} 2 . Let T e x 6x6 be the matrix taking v = (x 2 ,y 2 ,z 2 ,xy,xz,yz) T to 
(£i£[, i.2^21 ^3^3) Qii Q5, %) T ■ Then G = T T GT is the desired rank-3 Gram matrix 
for /, and all rank-3 Gram matrices arise in this way. Note that G does not 
depend on the choice of g 4 , q 5 , q 5 . □ 

Remark 5.11. Given /, finding a Steiner complex as input for the above algo- 
rithm is not a trivial task. But when a linear determinantal representation of / is 



M 



known, and thus a Cayley octad, one can use the criterion in Theorem 5.10 (3). 
Example 5.12. We consider the quartic / = det(M) defined by the matrix 

/52a; + 12y - 60z -26x - Qy + 30z 48z A8y \ 

-26x-Qy + 30z 2Qx + Qy - 30z -Qx+Qy-30z -A5x-27y-21z 
48z -6x + Qy-30z -9Qx A8x 

\ A8y -A5x-27y-21z A8x -A8x ) 

The complex curve Vc(/) is smooth and its set of real points Vr(/) is empty. The 
corresponding Cayley octad consists of four pairs of complex conjugates: 



6 1 



Here the 8x8 bitangent matrix OMO T = (6y) is defined over the field K = Q(i) 
of Gaussian rationals, and hence so are all 63 Gram matrices. According to the 



/ • 


—i 








-Q + Ai 


— 6 — 4i 


3 + 2i 


3-2i \ 


l + i 


l-i 








-4 + 4i 


-4-4i 


l-i 


7 + i 








i 


— i 


-3 + 2i 


-3 - 2% 


86 4 • 
39 13 


86 , 4 • 
39 ~ 13 


V 







l-i 


l-i 


l + i 


4 20 • 

39 39 


± + ^i 

39 ~ 39 / 
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lower right entry in Table [TJ precisely 15 of the Gram matrices are real, and hence 
these 15 Gram matrices have their entries in Q. For instance, the representation 



/ 



288 



(x 2 \ 


T 


/ 45500 


3102 


-9861 


5718 


-9246 


4956 \ 




( x2 \ 


y 2 




3102 


288 


-747 


882 


-18 


-144 




y 2 


z 2 




-9861 


-747 


3528 


-864 


-1170 


-504 




z 2 


xy 




5718 


882 


-864 


4440 


1104 


-2412 




xy 


xz 




-9246 


-18 


-1170 


1104 


11814 


-5058 




xz 






V 4956 


-144 


-504 


-2412 


-5058 


3582 J 




\v z ) 



is obtained by applying our algorithm for Theorem |5.1| to the Steiner complex 

S = {{^As}, {&15,&38}, {&18,&35}, {&24, h 7 } , {6 26) &47>, {&27, &46>}- 

The above Gram matrix has rank 3 and is positive semidefinite, so it translates 



into a representation ( 1.3 ) for / as the sum of three squares of quadrics over R. □ 



6. The Gram spectrahedron 

The Gram spectrahedron Gram(/) of a real ternary quartic / is the set of its 
positive semidefinite Gram matrices. This spectrahedron is the intersection of the 
cone of positive semidefinite 6x6-matrices with a 6-dimensional affine subspace. 
By Hilbert's result in [IB], Gram(/) is non-empty if and only if / is non-negative. 
In terms of coordinates on the 6-dimensional subspace given by a fixed quartic 

f(x, y, z) = C4002: 4 + c 3W x 3 y + c 301 x i z + c 2 2ox 2 y 2 + c 2 nx 2 yz H h c 00 4^ 4 , 

the Gram spectrahedron Gram(/) is the set of all positive semidefinite matrices 



(6.1] 



' C400 


Ai 


A 2 


2 C 310 


2 C 301 


A 4 \ 


Ai 


C040 


A3 


2^130 


A 5 


5C031 


A 2 


A 3 


C004 


A 6 


5C103 


^Coi3 


| c 310 


|ci30 


A 6 


C220 — 2Ai 


2 C 211 — A4 


|ci21 — A5 


| c 301 


As 


|ci03 


2 C 211 — A4 


C202 — 2A2 


\c\\2 — \% 


\ A 4 


5C031 


5C013 


3C121 — A5 


\c\\2 — A6 


CQ22 — 2A3 J 



where A G 



The main result of |24] is that a smooth positive quartic / has exactly eight 
inequivalent representations as a sum of three real squares, which had been con- 
jectured in |23]. These eight representations correspond to rank-3 positive semi- 
definite Gram matrices. We call these the vertices of rank 3 of Gram(/). In 
Section 5 we compute them using arithmetic over K. 

We define the Steiner graph of the Gram spectrahedron to be the graph on the 
eight vertices of rank 3 whose edges represent edges of the convex body Gram(/). 

Theorem 6.2. The Steiner graph of the Gram spectrahedron Gram(/) of a generic 
positive ternary quartic f is the disjoint union K4 U K4 of two complete graphs, 
and the relative interiors of these edges consist of rank-5 matrices. 

This theorem means that the eight rank-3 Gram matrices are divided into two 
groups of four, and, for G and G' in the same group, we have rank(G + G') < 5. 
The second sentence asserts that rank(G + G') = 5 holds for generic /. For the 
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proof it suffices to verify this for one specific /. This we have done, using exact 



arithmetic, for the quartic in Example 5.12 For instance, the rank-3 vertices 



-288^ (288^ 



/ 45500 3102 -9861 5718 -9246 4956 \ 

3102 288 -747 882 -18 -144 

-9861 -747 3528 -864 -1170 -504 

5718 882 -864 4440 1104 -2412 

-9246 -18 -1170 1104 11814 -5058 

y 4956 -144 -504 -2412 -5058 3582 J 



/ 45500 -2802 -6666 5718 -9246 132 \ 

-2802 288 -72 882 1206 -144 

-6666 -72 3528 -4878 -1170 -504 

5718 882 -4878 16248 5928 -3636 

-9246 1206 -1170 5928 5424 -1044 

y 132 -144 -504 -3636 -1044 2232 J 



both contain the vector (11355, -4241, 47584, 8325, 28530, 36706) T in their kernel, 
so that rank(G + G') < 5. But this vector spans the intersection of the kernels, 
hence rank(G + G') = 5, and every matrix on the edge has rank 5. 

We also know that there exist instances of smooth positive quartics where the 
rank along an edge drops to 4. One such example is the Fermat quartic, x A +y 4 +z 4: , 
which has two psd rank-3 Gram matrices whose sum has rank 4. We do not know 
whether the Gram spectrahedron Gram(/) has proper faces of dimension > 1 
other than the twelve edges in the Steiner graph K4UK4. In particular, we do not 
know whether the Steiner graph coincides with the graph of all edges of Gram(/). 



Proof of Theorem Fix a real symmetric linear determinantal representation 
M = xA + yB + zC of /. The existence of such M when / is positive was proved 
by Vinnikov [30l §0]. The Cayley octad {Oi, . . . , Og} determined by M consists 
of four pairs of complex conjugate points. Recall from Section [5] that a Steiner 
complex corresponds to either a subset / C {1, ... ,8} with |/| = 2 (type \J) or 
a partition I\I C of {1,. . . , 8} into two subsets of size 4 (type ||). We write Si 
for the Steiner complex given by / or I\I C and Gj for the corresponding Gram 



matrix. Theorem |6.2| follows from the more precise result in Theorem 6.3 which 



we shall prove further below. □ 

Theorem 6.3. Let f be positive with Vc(f) smooth and conjugation acting on the 
Cayley octad by Oi = O^a) for rr = (12) (34) (56) (78). The eight Steiner complexes 
corresponding to the vertices of rank 3 of the Gram spectrahedron Gram(/) are 

1357|2468 1368|2457 1458|2367 1467|2358 
1358|2467 1367|2458 1457|2368 1468|2357 

The Steiner graph K± U K4 is given by pairs of Steiner complexes in the same row. 



Our proof of Theorem 6.3 consists of two parts: (1) showing that the above 
Steiner complexes give the positive semidefinite Gram matrices and (2) showing 
how they form two copies of K 4 . We will begin by assuming (1) and proving (2): 



By Theorem 5.10 , for any two pairs of bitangents {£1, £[} and {£2^2} m a fixed 
Steiner complex S, there is a conic u in P 2 that passes through the eight contact 
points of these four bitangents with Vc(/)- In this manner, one associates with 
every Steiner complex S a set of (2) = 15 conies, denoted conics(5). 

Lemma 6.4. Let S and T be Steiner complexes with Gram matrices Gs and Gj-. 
J/conics(«S) fl conics(T) 7^ then rank(G lS + G T ) < 5. 
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Proof. Suppose S = {{£i,£[}, . . . , {A;^}}- Let Q be a quadratic matrix repre- 



sentation (5.5) such that the six points £i£' 1: . . . , £ 6 £' 6 G P(C[x, y, z\z) he on the 

we 



conic {X T QX : X G P 1 (C)}. By the construction in the proof of Theorem 5.1 
know that the projective plane in P(C[x, y, z]^) spanned by this conic is ker(G ? 5)" L . 

Consider two pairs {£i,£[}, {£ 2 ,£' 2 } from S and let u G conics(iS) be the unique 
conic passing through the eight contact points of these bitangents with the curve 
Vc(/)- By our choice of Q, we can find \/i G P 1 such that X T QX = £\£\ and 



^Q[i = £ 2 £ 2 . Equation (5.8) then shows that u = X^Qfi. From this we see that 
u G spanjQn, Qi 2 , Q22} = ker(G l5 )- L . Therefore, conics(iS) C ker(Gs) ± . 

If conics(iS) D conics(T) 7^ 0, then the two 3-planes ker(G lS )- L and ker(G 7 -) ± 
meet nontrivially. Since C[x, y, z] 2 has dimension 6, this implies that ker(Gs) and 
ker(GV) meet nontrivially. Hence rank(G,s + Gf) < 5. □ 

For example, conics(5i35s) and conies (51457) share the conic going through the 



contact points of 615, 626, &38, and 647. Lemma 6.4 then implies rank(Gi358 + 
G1457) < 5, as shown above for Example 5.12 with G = G1358 and G' = G1457. 



Using this approach, we only have to check that conics(iSf) D conics(5j) 7^ 
when / and J are in the same row of the table in Theorem |6.3| More precisely: 

Lemma 6.5. Let I and J be subsets of {1, ... , 8} of size four with I ^ J and 
I 7^ J c . Then conics(iS/) D conics(iSj) 7^ if and only if \I H J| = 2. 

Proof. Every syzygetic set of four bitangents £i,£ 2 ,£3, £4 determines a unique conic 
u passing through their eight contact points with Vc(/)- There are three ways to 
collect the four bitangents into two pairs, so u appears in conies (S) for exactly 
three Steiner complexes. Thus for two Steiner complexes Sj and Sj, we have 
conies (5/ )fl conies 7^ if and only if there are bitangents £\,£ 2 , £3, £4 such that 
{4,4}, {£3 A} e 5 7 and {4,4}, {4,4} € Sj. This translates into \I(lJ\ =2. □ 



To complete the proof of Theorem 6.3, it remains to show that the eight listed 



Steiner complexes give positive semidefinite Gram matrices. Recall that a Steiner 
complex Si is real if and only if / is fixed by the permutation rr coming from 
conjugation. As stated in Section |3j there are 15 real Steiner complexes, namely, 



(1) The eight complexes of type || listed in Theorem 6.3 

(2) Three more complexes of type ||, namely 1234|5678, 1256|3478, 1278|3456. 

(3) Four complexes of type \/, namely 12, 34, 56, 78. 

Since we know from |24j that exactly eight of these give positive semidefinite 
Gram matrices, it suffices to rule out the seven Steiner complexes in (2) and (3). 
Every Steiner complex Si gives rise to a system of contact conies Qi = {X T QiX, 



X G P 1 (C)}, where Qi is a symmetric 2x2-matrix as in (5.5), and a rank-3 Gram 



matrix Gi for /. The following proposition is a direct consequence of Remark 5.9 



Proposition 6.6. Let Si be a real Steiner complex. The Gram matrix Gi is 
positive semidefinite if and only if the system Qi does not contain any real conies. 

It follows that if Si is one of the three Steiner complexes in (2), then the Gram 
matrix Gj is not positive semidefinite, since the system Qi contains a product of 



20 



DANIEL PLAUMANN, BERND STURMFELS, AND CYNTHIA VINZANT 



two of the real bitangents 6 12 , 634, 656, bj$. Thus it remains to show that if / = ij 
with ij G {12, 34, 56, 78} as in (3), then the system Q^- contains a real conic. 

The symmetric linear determinantal representation M gives rise to the system 
{A T M adj A I A G P 3 (C)} of (azygetic) contact cubics (see H §6.3]). The main idea 
of the following is that multiplying a bitangent with a contact conic of / gives a 
contact cubic, and if both the bitangent and the cubic are real, then the conic 
must be real. The next lemma identifies products of bitangents and contact conies 
inside the system of contact cubics given by M. 

Lemma 6.7. For i ^ j we have ■ Q {j = {A T M adj A | A G span{0;, 0^}^}. 

Proof. After a change of coordinates, we can assume that Oi,Oj,Ok,Oi are the 
four unit vectors e\, ei, £3, This means that M = xA + yB + zC takes the form 



M 



(0 


bij 


bik 


bu\ 


bij 





bjk 


bji 


b ik 


bjk 





hi 


\bu 


b ji 


hi 


0/ 



Consider the three 3 x 3-minors complementary to the lower 2 x 2-block of M. They 
are e^M^e^, e^M^e^, ejM^e^. We check that all three are divisible by bij. 
Therefore 6" 1 • A T M ad jA with A G span{e3, e^} is a system of contact conies. Note 
that bikbjk = 6~ 1 eJ , M ad -'e4. Similarly, we can find the other six products of pairs 
of bitangents from the Steiner complex Sij, as illustrated by the following picture: 




Hence the system of contact conies Qi 



arises from division by b^ as asserted. 



□ 



Proof of Theorem 6.3 (and hence of Theorem 6.2). With all the various lemmas 
in place, only one tiny step is left to be done. Fix any of the four Steiner complexes 
ij of type \J in (3). Then the bitangent bij is real. Since M is real and Oi = Oj, we 
can pick a real point A G span{Oj, Oj} ± . Lemma 6.7 implies that that Qij contains 
the real conic b^f ■ A T M ad jA. Proposition 



6.6 



now completes the proof. 



□ 



Semidefinite programming over the Gram spectrahedron Gram(/) means find- 
ing the best sum of squares representation of a positive quartic /, where "best" 
refers to some criterion that can be expressed as a linear functional on Gram ma- 
trices. This optimization problem is of particular interest from the perspective 
of Tables 1 and 2 in |2Uj . because m = n = 6 is the smallest instance where the 
Pataki range of optimal ranks has size three. For the definition of Pataki range 
see also equation (5.16) in [25, §5]. The matrix rank of the exposed vertices of a 
generic 6-dimensional spectrahedron of 6x6-matrices can be either 3, 4 or 5. 

The Gram spectrahedra Gram(/) are not generic but they exhibit the generic 
behavior as far as the Pataki range is concerned. Namely, if we optimize a linear 
function over Gram(/) then the rank of the optimal matrix can be either 3, 4 or 
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5. We obtained the following numerical result for the distribution of these ranks 
by optimizing a random linear function over Gram(/) for randomly chosen /: 



Rank of optimal matrix 


3 


4 


5 


any 


Algebraic degree 


63 


38 


1 


102 


Probability 


2.01% 


95.44% 


2.55% 


100% 



TABLE 2. Statistics for semidefmite programming over Gram spectrahedra. 



The sampling in Table [2] was done in matlatj^J using the random matrix gen- 
erator. This distribution for the three possible ranks appears to be close to that 
of the generic case, as given in [20, Table 1]. The algebraic degree of the optimal 
solution, however, is much lower than in the generic situation of [2"0"l Table 2], 
where the three degrees are 112, 1400 and 32. For example, while the rank-3 lo- 
cus on the generic spectrahedron has 112 points over C, our Gram spectrahedron 
Gram(/) has only 63, one for each Steiner complex. 

The greatest surprise in Table [2] is the number 1 for the algebraic degree of 
the rank-5 solutions. This means that the optimal solution of a rational linear 
function over the Gram spectrahedron Gram(/) is Q-rational whenever it has 
rank 5. For a concrete example, consider the problem of maximizing the function 



159Ai - 9A 2 + 34A 3 + 73A 4 + 105A 5 + 86A 



<> 



over the Gram spectrahedron Gram(/) of the Fermat quarti c f = x +y +z . The 



optimal solution for this instance is the rank-5 Gram matrix (6.1 ) with coordinates 



-867799528369 -7785115393679 -2624916076477 1018287438360 2368982554265 562671279961 \ 



6890409751681 13780819503362 6890409751681 6890409751681 6890409751681 6890409751681 J 

The drop from 1400 to 38 for the algebraic degree of optimal Gram matrices of 
rank 4 is dramatic. It would be nice to understand the geometry behind this. We 
finally note that the algebraic degrees 63, 38, 1 in Table [2] were computed using 
Macaulay^Jby elimination from the KKT equations, as described in [25, §5]. 

7. The Variety of Cayley Octads 

The Cayley octads form a subvariety of codimension three in the space of eight 
labeled points in P 3 . A geometric study of this variety was undertaken by Dol- 
gachev and Ortland in [5; §IX.3], building on classical work of Coble [2j. This 
section complements their presentation with several explicit formulas we found 
useful for constructing examples and for performing symbolic computations. Be- 
sides convex algebraic geometry [T3l[T4*l 25J. our results have potential applications 
in number theory (e.g. arithmetic of del Pezzo varieties §V]) and integrable 
systems (e.g. 3-phase solutions to the Kadomtsev-Petviashvili equation [6J). In 



Theorem 7.5 we compute the discriminant of the quartics (1.2) and (1.3), and in 



Proposition 7.8 we discuss an application to nets of real quadrics in P 3 . 



2 

www . mathworks . com 



3 www . math . uiuc . edu/Macaulay2 
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We begin with the fact that a Cayley octad is determined by any seven of its 
points. Here is a rational formula for the eighth point in terms of the first seven. 

Proposition 7.1. Consider a general configuration C of seven points in P 3 , with 
coordinates (1:0:0:0), (0:1:0:0), (0:0:1:0), (0:0:0:1), (1:1:1:1), (a 6 :0 6 :^ 6 :5 6 ) and 
(a 7 :0 7 :^ 7 :b~ 7 ). The unique point (as'-Ps'-ls'-^s) ^ n which completes C to a Cayley 
octad is given by the following rational functions in the eight free parameters: 

0&17 - 0^7 - 1&07 + 76^7 + k07 - 5^ 7 

a 8 
0s 
7s 
5 8 

a 6 6 a 7 -f 7 - a 6 6 7 ^ 7 - a 6 j 6 a 7 7 + a 6 j 6 7 j 7 + 0e^a 7 7 - 0^%a 7 ^ 7 
Proof. This can be verified using linear algebra over the rational function field K = 
Q(«6, 06, 76, ^6, «7, 07, 77, fy) • We compute three linearly independent quadrics 
that vanish at the seven points and check that they also vanish at {a.%.0^%.8^). D 



^676^7^7 


~ /3e7677^7 - 06&60717 + 06^17^7 + 1^60717 ~ 
«677 - «6^7 - 76«7 + 76^7 + 0~6&7 ~ #677 


76^6 A^7 ' 




- "67677^7 - a 6 5 6 a 7 77 + a & 5^ 7 5 7 + ^5 & a 7 ^ 7 - 
a 6 7 - a 6 5 7 - 6 a 7 + 6 5 7 + 5 e a 7 - 5 6 7 


- 7 6 5 6 a 7 5 7 ' 


a 6 6 a 7 5 7 


- a 6 6 7 5 7 - a 6 5 6 a 7 7 + a 6 5 6 (3 7 5 7 + 6 5 6 a 7 7 
a 6 7 - a 6 77 - G a 7 + 0^ 7 + 7 6 a 7 - 7 6 /3 7 





The formula in Proposition |7. 1 1 parametrizes the semialgebraic set of real quar- 
tics that consist of four ovals. Indeed, if the parameters 0e, 76, Sq, «7, 07, 77, 87 



are real numbers then the corresponding quartic curve (1.2) has 28 real bitangents, 
so it falls into the first row of Table [TJ and all quartics with four ovals arise. Note 
that this row is the one relevant for applications to periodic water waves |6|. In 
practice we usually choose rational numbers for the parameters. This represents 



all curves whose 28 bitangents are rational, such as the Edge quartic (1.5), and it 
ensures that the ground field is K = Q for all computations in Sections 3-6. 

Yet the above formula has two disadvantages. First of all, it breaks the sym- 
metry among the eight points in the Cayley octad, and, secondly, it does not offer 
an arithmetically useful parametrization for the last two rows of Table [TJ Indeed, 
Vinnikov quartics and positive quartics are the lead actors in this paper, and we 



found ourselves unable to manipulate them properly using Proposition 7.1 For 



example, for a long time we failed to find a quartic with eight rank-3 Gram ma- 



trices over Q. Then we derived Proposition 7.2, and this led us to Example 5.12 
Let O be a configuration of eight points in general position in P 3 , represented 
by a 4x8-matrix. If O* is another such matrix whose row space equals the kernel 
of O then the configuration represented by O* is said to be Gale dual or associated 
to O. We refer to [SI IH] for the basics on Gale duality in the context of algebraic 
geometry. Both configurations O and O* are understood as equivalence classes 
modulo projective transformations of P 3 and relabeling of the eight points. We 
say that the configuration O is Gale self-dual if O and O* are equivalent in this 
sense. By a classical result due to Coble O is Gale self-dual if and only if O is 



a Cayley octad; see Theorem 3.10 



The variety of Cayley octads is defined by the equation = 0*. We now trans- 
late this equation into an algebraic form that is useful for computations. Let Pijki 
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denote the 4x4-minor of the 4x8-matrix that represents our configuration of 
eight points in P 3 . Consider the condition that O is mapped to a configuration pro- 
jectively equivalent to its Gale dual O* if we relabel the points by the permutation 
(18)(27)(36)(45). We express this condition using the Pliicker coordinates Piju- 

Proposition 7.2. Eight points in P 3 form a Cayley octad if and only if 



P1234P1256P3578P4678 
P1234P1267P3568P4578 
P1234P1457P2568P3678 
P1235P1267P3468P4578 
P1235P1367P2468P4578 
P1236P1347P2458P5678 
P1245P1267P3468P3578 
P1245P1356P2378P4678 
P1245P1367P2368P4578 
P1246P1357P2458P3678 



P5678P3478P1246P1235 , 
P5678P3458P1247P1236 j 
P5678P2368P1347P1245 , 
P4678P3458P1257P1236 , 
P4678P2458P1357P1236 j 
P4578P2568P1367P1234, 
P3678P3458P1257P1246 , 
P3678P2478P1456P1235 , 
P3678P2458P1457P1236 , 
P3578P2468P1367P1245 , 
and 



P1234P1257P3568P4678 
P1234P1356P2578P4678 
P1234P1467P2568P3578 
P1235P1347P2468P5678 
f>1235Pl467£>2468£>3578 
P1236P1456P2478P3578 
P1245P1346P2378P5678 
P1245P1357P2368P4678 
P1246P1357P2368P4578 
P1247P1357P2368P4568 
P1346P1357P2458P2678 



P5678P3468P1247P1235 , 
P5678P2478P1346P1235 j 
P5678P2358P1347P1246 ; 
P4678P2568P1357P1234 , 
P4678P2358P1357P1246 j 
P4578P2378P1356P1246 ; 
P3678P2578P1456P1234 , 
P3678P2468P1457P1235 j 
P3578P2468P1457P1236 ; 
P3568P2468P1457P1237 , 
P2578P2468P1367P1345 • 



Before discussing the proof of this theorem, we first explain why the shape of 
the above equations is plausible. Consider the condition for six points (xj : yi : Zi) 
in P 2 to be self-dual, in the sense above. This condition means that the six points 
lie on a conic, and we write this algebraically in terms of Pliicker coordinates as 



(7.3) det 



(x\ 



X2 
4 

4 

We 



v\ 
vi 
vl 
vl 
vl 

vl 



,2 
z l 

"2 

~3 

z A 

■y'2 



xiyi 

X2V2 

X5V5 
X6V6 



X\Z\ 
X2Z2 
X3Z3 
X4Z4 

X5Z5 
X 6 Z 6 



yizi\ 

V2Z2 
V3Z3 
ViZi 

V5Z5 
V&ZqJ 



Pl23£>145£>246£>356 ~ ^124^135^236^456 • 



This formula appears in [5j Ex.4, p. 37] and we adapt the derivation given there. 



Sketch of proof for Proposition 7.2. The cross ratio (]?i234Pi256)/(Pi235Pi246) is in- 
variant under projective transformations. The permutation (18) (27) (36) (45) of 
the points transforms that cross ratio into (ps765P8743)/(P8764P8753)- The condition 
= 0* implies that these two cross ratios are equal. By clearing denominators, 
the equality of cross ratios translates into the first of the 21 equations listed above: 



Pl234£>1256P3578P4678 



^5678^3478^1246^1235 



The other 20 equations are found by the same argument for cross ratios. By 
incorporating the quadratic Pliicker relations among the Pijki, we check that our 
list of 21 cross ratio identities is complete, in the sense that it ensures = 0*. □ 

Dolgachev and Ortland |5} page 176] present the conditions under which a 
regular Cayley octad O can degenerate. Their analysis exhibits 64 = 28 + 35 + 1 
boundary divisors in the compactified space of Cayley octads. These are as follows: 

(1) Two points of O can come together. This gives 28 = (2) boundary divisors, 
e.g., points 1 and 2 come together if and only if p\uj = for 3<z<j<8. 
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(2) Four points of can become coplanar. The equations in Proposition 7.2 
then ensure that the other four points become coplanar as well. So, in 
total there are 35 = |(J boundary divisors such as {^1234 = P5678 — 0}. 

(3) The eight points of O can lie on a twisted cubic curve, which is the inter- 
section of the three quadrics. The condition for seven points in P 3 to lie on 
a twisted cubic curve has codimension 2. H.S. White [2U| eqn.(2)] writes 



this condition by adding an index to (7.3). This gives 7-15-3 equations like 

(7-4) P1237P1457P2467P3567 - P1247P1357P2367P4567 = 0. 

Applying the symmetric group Sg to the indices, we obtain equations for 
the codimension 4 locus of octads that lie on a twisted cubic curve. This 
locus is a divisor in the compactified space of Cayley octads, as in 0, 



5IX.3]. Equivalently, the equations (7.4) imply those in Proposition 7.2 



We now shift gears and examine the three types of boundary divisors from the 



perspective of the desirable representations (1.2) and (1.3) of a ternary quartic /. 
In other words, we wish to identify the conditions, expressed algebraically in terms 
of these two representations, for the quartic curve Vc(/) to become singular. 
Recall that the discriminant A of / is a homogeneous polynomial of degree 27, 



featured explicitly in [271 Proposition 6.5], in the 15 coefficients Cijk of (1.1). If 



we take / in the representation (1.2) then each coefficient c^k is replaced by a 
polynomial of degree 4 in the 30 = 10 + 10 + 10 entries of the symmetric matrices 
A, B and C. The result of performing this substitution in the discriminant A(cijk) 
is denoted A(A, B,C). This is a homogeneous polynomial of degree 108 in 30 
unknowns. We call A(A,B,C) the Vinnikov discriminant of a ternary quartic. 



Similarly, if we take / in the representation (1.3) then each coefficient Cy& is 
replaced by a polynomial of degree 2 in the 18 = 6+6+6 coefficients of the quadrics 
91, q 2 and 93. The result of performing this substitution in the discriminant A(cjjfe) 
is denoted A(qx, 92,93)- This is a homogeneous polynomial of degree 54 in 18 
unknowns. We call A(qi, q2, 93) the Hilbert discriminant of a ternary quartic. 

Theorem 7.5. The irreducible factorization of the Vinnikov discriminant equals 

(7.6) A(A,B,C) = M(A,B,C) -P{A,B,C) 2 , 

where P has degree 30 and corresponds to the boundary divisor (2), while M has 
degree 48, and this is the mixed discriminant corresponding to both (1) and (3). 
The irreducible factorization of the Hilbert discriminant equals 

(7.7) A{q u q 2 ,qz) = Q(q u q 2 , q 3 ) • R(?i, 92, 93) 2 , 

where Q has degree 30 and the degree 12 factor R is the resultant of q±, q 2 and 93. 

This theorem is proved by a computation, the details of which we omit here. 
It has been pointed out to us by Igor Dolgachev and Giorgio Ottaviani that the 



factorization (7.6) was already known to Salmon [26J, who refers to ~M.(A,B,C) 



as the tact invariant. See also |TXJ| Section 10] for a modern treatment. 



We discuss the geometric meaning of the factors in (7.6) and (7.7). The poly- 



nomials M, P, Q, and R are absolutely irreducible: they do not factor over C. 
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The polynomial P represents the condition that the span of A, B and C in the 
space of 4x4-symmetric matrices contains a rank-2 matrix. Note that the variety 
of such rank-2 matrices has codimension 3 and degree 10. The Chow form of that 
variety is precisely our polynomial P, which explains why P has degree 3- 10 = 30. 

Non- vanishing of the mixed discriminant M is the condition for the intersection 
of three quadrics in P 3 to be zero-dimensional and smooth. A general formula for 
the degree of such discriminants appears in [T9j Theorem 3.1]. It implies that M 
is tri-homogeneous of degree (16,16,16) in the entries of (A,B,C), so the total 
degree of M is 48. Note that vanishing of M represents not just condition (1) but 
it also subsumes condition (3) that the quadrics intersect in a twisted cubic curve. 

The resultant R of three ternary quadrics q 2 , q^) is tri-homogeneous of de- 
gree (4,4,4) since two quadrics meet in 4 points in P 2 . Thus R has total degree 
12. The extraneous factor Q of degree 30 expresses the condition that, at some 
point in P 2 , the vector (qi, q%, q?) is non-zero and lies in the kernel of its Jacobian. 

We close this paper by reinterpreting Table [I] as a tool to study linear spaces of 
symmetric 4x4 matrices. Two matrices A and B determine a pencil of quadrics in 
P 3 , and three matrices A, B, C determine a net of quadrics in P 3 . We now consider 
these pencils and nets over the field K of real numbers. A classical fact, proved 
by Calabi in |T], states that a pencil of quadrics either has a common point or 
contains a positive definite quadric. This fact is the foundation for an optimization 
technique known in engineering as the S-procedure. The same dichotomy is false 
for nets of quadrics [H, §4], and for quadrics in P 3 it fails in two interesting ways. 

Theorem 7.8. Let Af be a real net of homogeneous quadrics in four unknowns 
with A(jV) 7^ 0. Then precisely one of the following four cases holds: 

(a) The quadrics in Af have a common point in P 3 (IR). 

(b) The net Af is definite, i.e. it contains a positive definite quadric. 

(c) There is a definite net Af' with det(7V') = det(jV) ; but Af is nondefinite. 

(d) The net Af contains no singular quadric. 

Proof. For a real net of quadrics, Af = $L{A,B,C}, the Vinnikov discriminant 
A(A,B,C) in (7.6) is independent (up to scaling) of the basis {A, B,C}, and 
thus can be denoted A(Af). If A(Af) is non-zero, the polynomial det(AT) = 
det(xA + yB + zC) defines a smooth curve, which depends on the choice of basis 
{A, B,C} only up to projective change of coordinates in [x : y : z\. This real 
quartic falls into precisely one of the six classes in Table [T] The first four classes 
correspond to our case (a). The fifth class corresponds to our cases (b) and (c) 
by the Helton- Vinnikov Theorem [13J. As a Vinnikov quartic has definite and 
non-definite real determinantal representations, both (b) and (c) do occur [30J. 
For an example, see [2"T| Ex. 5.2]. The last class corresponds to our case (d). □ 

Given a net of quadrics Af = M{A, B,C}, one may wish to know whether there 
is a common intersection point in real projective 3-space P 3 (M), and, if not, one 
seeks the certificates promised in parts (b)-(d) of Theorem |7.8| Our algorithms 
in Sections 3, 4 and 5 furnish a practical method for identifying cases (b) and (d). 
The difference between (b) and (c) is more subtle and is discussed in detail in |2"T] . 
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